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Thermodynamics of layered Ising magnet with the infinite-range ferromagnetic intralayer interac- 
tion and random exchange between nearest layers is considered. The case of zero average interlayer 
exchange is studied in detail. The inequilibrium thermodynamic potential is obtained near transi- 
tion point and existence of numerous metastable states is shown. The thermodynamic properties of 
crystal in the simple periodic states are described. It is established that the evolution of the equi- 
librium magnetic state with growth of random exchange proceed over infinite series of first order 
transitions at which the number of magnetized layers and distance between them jump as well as 
homogeneous magnetization and susceptibility. 
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(— I The complete or partial disorder of magnetic structure has been found in many layered magnets and to describe 

it the term 'spin glass' has often been usedi'2iiiiSi&. The other spin glass features were also observed in layered 
I , magnets such as the difference of 'field-cooled' and 'zero-field-cooled' parameters^afe, smeared peak and low-frequency 
, ^ ■ dispersion of magnetic susceptibilitji^i^. The glassy properties do not exist only in the crystalline solid solutions, in 
O which they obviousely result from magnetic exchange fluctuations caused by ion substitutions. Some crystals allow 
significant deviation from stechiometry, such as SrCoOs-s^, in others, the disorder is their intrinsic property as in 
LiNi02 - here Li and Ni ions can easily exchange positions as they have almost equal sizes ^. Also the nominally 
pure crystals has always some impurities and defects. 

The simplest variant of disorder is the intercalation of impurities into the interlayer space, which does not modify 
fH , the intralayer magnetic interactions. Similar model explains the glassy features of LiNi02 - small concentration of 
Q ' Li ions in the ferromagnetic layers containing Ni does not essentially affect them, while the presence of Ni ions in Li 
O , layers results in appearance of local ferromagnetic interlayer exchange on the background of small antiferromagnetic 
onc''^. In this case one may hope to get comparatively simple description of disorder effects, especially when radius of 
^ magnetic intralayer interaction is much greater than a lattice parameter. Then intralayer ordering can be considered 
^ . in the mean-field approximation everywhere except a narrow vicinity of transition point. Here the interlayer exchange 
lO ' fluctuations are averaged over intralayer interaction range, which can be put infinite in this mean-field region. Hence, 
we get a quasi one-dimensional situation with some average interlayer interaction''. Then theoretical description 
■ becomes simpler than in case of 3d disorder and one may try to get some analytical results on the thermodynamics of 
such random layered magnet, preserving some specific features of spin glasses, particularly, the presence of multiple 
metastable states. The origin and properties of such states is a central issue of the theory of random magnets, so the 
' study of such model could give not only approximate description of thermodynamics of real crystals but also shed 
. some light on the mechanism of emergence of multiple metastable states responsible for the inergodic behavior of spin 
d ' glasses. 

^ , Here we consider the simplest model of Ising magnet having infinite-range ferromagnetic intralayer interaction and 
I random exchange between spins in neighboring layers. It allows to obtain rather easily an effective thermodynamic 
' ^ , potential for the layers' magnetizations near transition point. Therefore, we can make sufficiently complete study of 
the origin and the properties of various metastable states for the case of zero average interlayer exchange. Further we 
establish the existence of a number of equilibrium spin-glass phases for different values of random exchange fluctuations 
and describe phase transitions between them. 
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I. MODEL AND FORMALISM 

The Hamiltonian of the model is 

J N / M \^MN M N 

^ '2 A I ^ ^ ( ^ ^ ^ra,n 1 ^ ^ ^ ^ Jm,7iSm,nSm,n+l ^ ^ ^ ^ Srn.nHn (1) 

n—1 \m— 1 / ni—1 n—1 m—1 n—1 

Here Sm.n = ±1 are Ising spins, index m numerates the lattice sites in the layer from 1 to M, n = 1,2,..., TV is the 
layer number, J > is ferromagnetic intralayer exchange, Jm,n describes random interlayer interaction and Hn is 
magnetic fleld, which is constant inside a layer. We assume all Jm,n to have the same distribution function. 
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For the Hamiltonian in Eq. one can get the averaged over Jm.n value of the thermodynamic potential (per site 
in a layer) 

MPF = - (InTr exp {-f3n))j . 

Here fi = 1/T is inverse temperature, Tr means the sum over all spin configurations and the angle parenthesis with 
index J below designate the average over random Jm,n- 

The Hubbard-Stratonovich transformation of the ferromagnetic term in Eq. gives 



M/3F = - ( In <^ dtcxp -Ml3F (f 



(2) 



Af ^ M f N 

n—1 m—1 \ n—1 



m,n+l ~t~ {hji H" /n) Sjn^n 



(3) 



where /i„ = /3i7„ , K ^ pj and K,n,n = PJm,n 

As the logarithms in the sum in Eq. (PJ are independent random quantities the averaging of Eq. Q over random 
interlayer exchange reduces to the averaging of this large sum when M —^ oo. Thus we get 



Mf3F ^ -liii d{ exp (f )] 



N 



In 



Tr 

n=l ) I J 



(4) 
(5) 



In Eq. ^ index m is omitted at spins and random exchange integral because the contributions from all sites 
in the layers are identical. Obtaining the inequilibrium thermodynamic potential in Eq. |(SJ), one can describe all 
thermodynamics of the model. Indeed, at M — > oo we have 



F « mm [F (f )] - TSconf, 



(6) 



where Sconf the is configurational entropy defined via the logarithm of the number W of points f at which the 
minimum of F (f ) is attained 



For the layer magnetization 



M 



m—1 



dl3F 
dhn ' 



we get from Eqs. (01 



/ dl3F (f ) \ 



d(3F (f ) 



and for the correlator defining the inhomogeneous susceptipility 

M 



(7) 



d'^PF 



dhndhn' 



we have 



Cn,„, = MR- 



(/ra/n')/ {fn)f {fn') f 



^d'^lSF (f) 



didi 



f=(f) 



(8) 
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The second term in Eq. ((SJ contains under the logarithm the partition function of Ising chain with random exchange 
in inhomogeneous field. It can be represented as expansion in powers of /i„ = tanh (/i„ + /„) and one can express the 
potential in Eq. © as a function of these variables 



pF ill) =-{N-l) (\n (2 coshi^) + [(tanh-i(Ai„) - h^f /K + In (l - /i') 

/ / lN/2] 2k \\ 

-y'^i^^Y. g„,,„,,...,„,,J|a^„, ) , (9) 



k—1 ni<n2<...<n2fc 




To perform the averaging in the last term of Eq. (j^ is generally a hard task, it is much simpler to study the vicinity 
of transition point and the case of small or zero external fields. Then we can leave in Eq. only the lowest powers 
in its expansion in small /x„ and Further we consider just this case. 



II. THERMODYNAMICS NEAR TRANSITION POINT 

Up to the forth order in /i„ and for |/i„| < we get from Eq. @ 

PF = - (iV - 1) (in (2 coshi^) ) _ + ^ ^ + i ^ - gW„,) + fJ- - /4 + 

n,n' ^ ' n 

1 2 (2) 2 , 2 (2) (1) , (1) (2) 2\ , o Y^ (1) (2) (1) 

2Z^M^5^->j+ 1^ \y'i9i-jH9)-ktJ'k + ^J'^9i-J^J'j9)-k^J'k) +^ 2^ tJ'i9i-ji^i9)lky'kglliiJ^i [iQ) 

i<j i<j<k i<j<k<l 

where gn'^ — v''^', gn'' = u'"', v = (tanh/?^ ^, u = (tanh^ k'^ 

As Hn ~ fn/K, the equilibrium (corresponding to the potential minimum) values of these variables give the layers' 
magnetizations, cf. Eq. (UJ. Apparently, in zero field the potential in Eq. H10|l describes the second order phase 
transition at which spontaneous layers' magnetizations order ferromagnetically when w > or antiferromagnetically 



(1) 

gti 



= 1. 



when V < 0. The transition point in both cases is defined by the condition K ^ 

n— — i 

In general, the essentially nonlocal interaction of the forth order in Eq. (|1U|) (its radius goes to infinity when 
w — *■ 1) allows to suppose that along with the minimum of the type /x„ = /x or /i„ = (—1) the other minima could 
exist corresponding to some inhomogeneous ordering. Most probably, these minima would have larger values of the 
potential, thus defining the metastable states of a crystal. 

Yet, in the presence of a field (especially inhomogeneous one) these states could become stable, i. e. there would 
be first order transitions from ferro- or antiferromagnetic state into various inhomogeneous ones depending on the 
strength and spatial variations the field applied. The dynamic phenomena in slow varying external fields would 
also have pecularities caused by the presence of inhomogeneous metastable states so their study in the ferro- or 
antiferromagnetic cases are very important for the description of the observed glassy properties of layered magnets. 
The interesting features are also revealed in the numeric study of temperature dependence of the magnetization in the 
case of random ferromagnetic excange^ - it shows visible oscillations growing while temperature is decreasing below 
the transition point. Nature of this phenomenon is not clear now. 

However, the case of pure glass ordering, realized when the distribution function of interlayer exchange is symmetric, 
seems to be the most interesting one. Then v — and the potential in Eq. (|10|) at /i„ = becomes 

PF in) = - (iV - 1) (in (2 coshi^) )^ + ^A^' - ^ E ^" + i E (11) 

n n,n' 

where r = — 1 = (T — J) / J, r <C 1. In spite of apparent simplicity of Eq. I|1HI the determination of spontaneous 
magnetizations /i„ appearing at r < is rather complicated. The equation for the F (fi) extrema 

/^nf--^/^n+E"'"""A^M -0, (12) 
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has in general case a number of solutions, among which one should choose the stable ones (i. e. corresponding to 
local minima) having the positively definite matrix of second order derivatives (Hessian) 

= W (^r-2Ai?, + ^«l"-Vfj +2/i„«l"-"V„.. (13) 

The form of stable solutions of Eq. (|12|l depends essentially on the u value. Indeed, in the trivial case m = they are 
of the form /i„ — ±/x, i. e. all layers have the same magnetizations' magnitudes while their orientations are arbitrary, 
so there are 2^ solutions. This degeneracy of the potential minima is the consequence of its symmetry with respect 
to the sign change of every /x„. In the other extreme when random interactions are infinite, m — > 1, the stable states 
have only one (arbitrary) layer with nonzero magnetization, thus there are 2N solutions (with different fin ^ signs). 
Our task is to find out which states are realized at the intermediate u values and how the evolution of the magnetic 
state of crystal proceeds when disorder grows. 

It follows from Eq. ()12|l that general form of its solution is 



= i^nVFk^ (14) 
The variables can take two values, and 1, i?„ = {0, 1}. Choosing the definite vector -d, one can find q from Eq. 

m 

q = D-'d, (15) 

Dn,n' =^nu\"-'''\^n' ~^Sn,n'. (16) 

Let us note that variable |T|q has the meaning of the (inhomogeneous) Edwards- Anderson order parameter^. Ac- 
cording to Eq. 1)14(1 the choice of ■& for a given u is limited by the condition g„ > 0. Also the positive definiteness is 
necessary for Hessian, Eq. ()13|l . or, equivalently, for the correlator, cf. Eq. which has the following form on the 
solutions considered 



V / n.n' 



(17) 



For this the validity of the following conditions are sufficient 

En>0 at i?„ = 0, (19) 

^nDn^n'K' > 0. (20) 

The last inequality designates symbolically the positive defineteness of the matrix D on the subspace where = 1. 
The corresponding values of the equlibrium thermodynamic potential are 

2 

f3F = -{N ^1) (in (2 cosh if) ^^^ + _ g^^^^ ^ _ (^y - 1) (|ln (2coshi:) - ^i^ZJ-^i? - M^^ ^ In 2. 



III. PERIODIC STATES 



To find among 2 vectors those obeying the above conditions at a given u is rather difficult task. So we consider 
first only periodic states having some period L. Introducing the layers' numeration of the form n — rL -\- s with 
integer r and s=l,2,...,Lwe will study the solutions such that "drL+s = "^s- Apparently, the i-dimensional vectors 
z?s , which differ by a cyclic permutations of components, describe the same solution. Then it is convenient to consider 
"ds as defined on the equidistance points of a circle to show visually that the rotation of -dg on a circle do not change 
the state. 

Let us also note that the signs of /i„ in the states considered should not be necessarily periodic as they do not enter 
the stability conditions and the equilibrium potential value. So further results refer to all 21^/^1^ states differing by 

L 

layers' magnetization orientations. Here P = ^ i?s is the number of magnetized layers in a period and the quadrangle 

s=l 

parentheses denote the integer part of a number. 
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For such periodic states we have 
As the matrix D depends only on r — r' it can be diagonalized over this indexes using the unitary transformation 

1 /2 

fife^r = (L/N) ' expikr in which we assume that k belong to the first Brillouin zone, \k\ < tt. The result is 



1 u~^e 



2. 



(21) 



Now the stability condition for the matrix D, Eq. H2U|I . means that Dg^s' (k) must be positively defined on subspace 
where i?^ = 1 at all |fc| < tt. Periodic solutions for q can be expressed via Ds.s' (k) at A; = 0, 



QrL+s — Qs — 



D{oy 



the same is true for Eg 

ErL+s ^ Es = Ds,s' (0)(?s'l^s' - 1 

Diagonalization of the correlator over r — r' gives 



and for the equilibrium potential we have 



D{ky 



PF/N K - (\n{2coshkjj ^- 



-i?L»(0) ^i9 + M"ipin2 



(22) 



(23) 



(24) 



(25) 



The last term in Eq. H25|) proportional to P is the contribution of the configurational entropy appearing due to 
degeneracy of the potential minima as there are 2[^/^l-^ states with equal potentials differing by the signs of layers' 
magnetizations. It should be taken into account when P is of order M. As P < N, this could only happen ii N M, 
i. e. in thin whiskers with diameter of order \/M much smaller than the length A''. Yet we should note that this 
degeneracy are present only in the idealized situation when an external field is exactly zero while in real experiments 
there always is some small inhomogeneous field breaking this degeneracy. Further we assume the presence of such 
infinitesimal field and drop the contribution from the configurational entropy irrespective of the sample geometry. 



IV. THERMODYNAMICS AND STABILITY OF PERIODIC STATES 

Let us consider the simplest periodic states for which analytical results can be obtained. Thus, for the states with 
P — 1 {-Ol — \,'ds = 0, s = 1, L — 1) the matrix D in Eq. H21() becomes just the function of fc. 



and we have 



^ 11 + Au^ cos k - 5u'^^ 

^ ' ~ 3 1 - 2ui cos /c + 1*2^ 



qL = q = D{0) ^ " 



E.. 



u + u 



L-s 



1 + ' 
l,...,i- 1 



(26) 

(27) 
(28) 



Substitution of these expressions into Eq. 124|l gives the formulae for the correlator which allows to find, in 
particular, the expression for the homogeneous susceptibility. 



PX = (0) = ( 1 + 2 Xi / (2 |r| L) 

s,s' \ s=l / 



(29) 
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It follows from Eqs. f??. l20(l that the stability region of these states is defined by the conditions E^l/2\ > 0, D (tt) > 0. 
The first of them defines the lower bound for u values and the last one defines the upper bound. For the even L we 
have 

< (30) 

For small odd L the lower bound is slightly higher, but while L grows it soon comes to 1/25. 
Introducing the variable 

Riu) = -^^, (31) 
mu 

having a meaning of effective radius of interlayer exchange, the stability condition in Eq. I|3U|I can be represented in 
the form 

1<L/R<2. (32) 

This representation demonstrates explicitly that period L (the distance between magnetized layers) cannot significantly 
differ from the radius of interlayer exchange. The equilibrium potential of the states considered is 



PF/N (In [2 cosh Kj j^- r^q/AL. (33) 

Now we turn to the states with P = 2. Let -di^ — ^i^r — 1, while the rest i9s — and assume for definiteness 
L' < L/2. Then 

3(1-^-^) 

^^ = '1^'-^= , + 3,L^ + ^,L-L' + (34) 
PF/N ^ - (in (2 cosh K^"^ ^-T^q'/2L (35) 

For these states to be stable the conditions i?[(L+L')/2] > and Dll (0) — Dll' (0) > must be fulfilled or (for L + L' 
even) 



(^jL-L')/2 ^ „(L+L')/2) _ 3 (^^L' ^ ^L-L'^ 



- 5u^ - 1 > 0, 



1 + 5m^ - 3 (u^' + u^^^j > (36) 

Summing these inequalities we get the condition u^^'" ^^)/2 < 1 so they both can be satisfied only when L' > L/3 or 

l<iy <2 (37) 

where v = [L — L')/ L' is the ratio of the maximal distance between magnetized layers and the minimal one. 
The stability conditions in Eq. H36|l can be presented in the form similar to that of Eq. H32() 

A' [v] < L'/R < B' [v) . (38) 

Numerical analysis of Eq. ()36f) gives the following approximate expressions 

A' {v) « 0.77 + 0.23/1/, B' {v) w -0.23 + 2.2i/v 

For V — 2 the stability region reduces to a point while dX v = 1 {L' = L/2) it coincides with Eq. H32I) because the 
corresponding states are identical to those having P = 1 and Li = L/2. It is easy to see that these states with P = 1 
have the stability region (L/4 < i? < L/2) which includes the whole interval of Eq. L'/B' (ly) < R < L'/A' (ly). 
They also have the lower potential values than the P — 2 ones. The last statement is the consequence of the maximal 
value of g', Eq. H34|l . at L' — L/2. Thus the states with P = 2 are always metastable. 

Analytical study of states with P > 2 is very cumbersome, yet their main properties could be established using 
numerical methods and qualitative considerations following from the exact results for P = 1, 2. Thus one can deduce 
from the preceding results that at all P stable states should have the distances between magnetized layers of the order 
R. Introducing the minimal and the maximal values of such distances (and assuming the points representing layers 
to lie on a circle) 

imin = min(|s - s'l ,L - |s - s'l) , Lmax = max (|s - s'| , L - |s - s'|) for 'ds'ds' = l 
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FIG. 1: . The functions A{u) (quadrangles) and B{i') (circles) in Eq. 14111 obtained in the numerical studies of stable states 
with P > 2. 



one can get the upper bound for the parameter v = Lmax/-^min of the stable states. 

Indeed, one can see that the stabihty condition for the matrix D(0) in Eq. (|21|l imphes that its maximal nondiagonal 
element (on the subspace where ??s = 1) must be smaller than the diagonal ones. Considering states with P ^ 1 
(and, hence, L » R) we get from this the approximate necessary condition 

while from the condition Ei^^_^^/2 > 0, cLSSuming LmaK ~ R 0-iid ~ Qrnax — 3 ill Ec[. (|23|) we have 8,pproxiiiicLtely 

6yin,ax/2 > I (-40) 

From Eqs. H39I I40() it follows that only those states are stable in finite R interval which have 



« 3.262. 

In 3 



The largest f value found in numerical studies of states with 5 < L < 30 and 3 < P < 15 is 3. This corroborates the 
existence of upper bound i^c ~ 3 for stable states. Numerical results show also that for P > 2 there are the necessary 
stability conditions similar to those of Eq. H38|) , 



A (ly) < L^in/R <B{iy). 



(41) 



The functions A (i/) and B (v) obtained in the numerical studies are shown in Fig^ Note that the region between 
A (ly) and B (ly) includes the region defined by Eq. H38(l and that Eq. H41|) coincide with Eq. (|32|l at ^ 1 as in this 
case all states become equivalent to those with P — 1. From Eq. (|41|l it follows that interval of R values in which 
stable states with P > 1 and 1 < j/ < f c « 3 exist, 

imin/B (ly) < R< LminM i") > 

belongs entirely to the stability region of states with P = 1, Eq. H32(l . and some Li if Li obeys the condition 

Lnun/A (iy) <Li< 2Ln,in/B (z/) . 

Apparently such Li do exist and owing to A (v) < 1 and i^B (v) > 2 their values are confined between Lmin and £max, 

^min ^ Li < -/^max- 
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There are also states with Li in this interval and in its vicinity which cover partially the stability region of Eq. 1)41(1 . 

The mentioned above numerical studies show also that states with P > 1 have always the higher thermodynamic 
potential than the coexisting with them P = 1 states . The necessity of this can be shown analytically for the states 
with V close to 1. Let us consider the states with P > 1 and L = PLi having equidistant magnetized layers, i. e. 
equivalent to the states with P — 1 and Li — L/P. We can find the change of their potential, Eq. (|25() . under small 
variations of the positions of magnetized layers, Lis Lis + Srg. The result is 

PSF/N = {Tq\nu/2f \Srk\^ \d (0) ^ D [k) - D {ky^ \D' {k)f 

k 

Here 5rk is Fourier-transform of 6rs, k = 2ttI/P, I = 1,2, P, D (fc) is given by Eq. (|26|) with L = Li and 



2iui sin k _ ^ 

L> (k) = ; TT, Ul = u 

^ ^ 1 -2ui COS k + uj' 



In the stability region (ii/2 < R < Li), where D (k) > 0, the expression under the sum sign is positive for all k so 
6F > and states with P > 1 and h' close to 1 have always the higher thermodynamic potential than the coexisting 
with them P = 1 states. 

We cannot give a rigorous proof of the metastability of states with P > 1 at all < i/c- However, the numerical 
results allow to suggest that all states having different distances between magnetized layers have higher potential then 
that of the states with P = 1 in which only one such distance (L) exists. The last means that nonperiodic states will 
also be metastable so equilibrium thermodynamic properties of the model are solely determined by the states with 
P = 1. Further we consider the equilibrium thermodynamics and phase transition in the model in the framework of 
this supposition. 



V. PHASE TRANSITIONS AT DISORDER STRENGTH CHANGES 

The stability condition of P = 1 states, Eq. (|32|l . represented as R < L < 2R shows that there are several stable 
states at all P > 1. Hence one should find among them the states with the lowest potential, Eq. H33|) . to describe the 
equilibrium properties of a crystal. Considering the potential as continuous function of L we find that it has minimum 
at 

Lc^\R, Aw 1.471 

Thus the state with the lowest potential has L closest to Lc- At 

Pl = A-1(L + 1/2), L>1. 

there are two nearest and equidistant from Lc periods (L and L + 1), so at that the jumps of the distance between 
magnetized layers occur manifesting the first order phase transitions. At P = 1 there is also phase transition between 
L = 1 state (where all layers are magnetized) and L = 2 state, this point being the stability boundary of both these 
phases. 

Thus the evolution of crystal state under the growth of disorder fluctuations (increasing of u or P) proceed via 
infinite series of phase transitions accompanied by the jumps of the distance between magnetized layers as well as 
thermodynamic parameters of a crystal. To get this equilibrium parameters one should substitute in Eqs. (|27l 1291 . 
) the values of L corresponding to a state with the lowest potential 

oo 

Leq{R) = I? (1.5 - P) [1 + I? (P - 1)] + ?9 (P - 1.5) ^ Li? [l - 4 (AP - L) 



L=2 



This gives 



1 _ 5-L.,(fl)/i? 

Qeq (P) = 3- 



1 + 5l-L.,(fl)/i?' 



L,,(_R)-1 

I3xeq (P) = I 1 + 2^ (P - 1.5) J2 I / (2 \r\ ie, (i?)) 
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FIG. 2: . The dependences of geq (quadrangles), \r\l3xeq (circles) and 5C (triangles) on the effective interaction radius R. 
The jumps of first two parameters correspond to the phase transitions between the states with different distances between 
magnetized layers. 



Let us remind that Qeq (R) defines the spontaneous magnetization of layers nieq (R) — \/\t\ qeq (i?). The behavior of 
Qeq (R) and the equilibrium homogeneous susceptibility is shown in Fig[21 Differentiating the equilibrium potential 
over the temperature we get the magnetic contribution to the heat capacity appearing at r < (T < J) 

5C = qeq (R) /2Leq (R) . 

It has no jumps at the transitions being proportional to the equilibrium potential, see Fig[21 

It must be noted that these thermodynamic parameters are the same for all 2^^/^! states differing by the layers' 
magnetizations orientations while the total homogeneous magnetization is apparently different. It is determined by 
the infinitesimal field existing in the experiment and choosing the state in which > for all n. Then the 

overwhelming majority of these randomly chosen states would have zero or almost zero total magnetization. Indeed, 
the number of states with magnetization m is 

Nm = ( - m) /2 ) ^ exp [(1 + m) In (1 + m) + (1 - rn) In (1 - m)] 

It has narrow peak at m = with the width of order [N/ L] ' . 

The results obtained so far can be slightly changed if potential expansion in fi will include higher order terms. 
Primarily, it concerns the vicinity of i? = 1 point being the stability boundary of neighboring states. Here both 
phases becomes unstable with respect to the transition into L = 2 phase of the type [^1,^2) which is absent in 
the forth-order approximation considered. Considering the higher order terms could result either in appearance of 
narrow (of order r) interval of existence of such phase with second order transition points at its boundaries, or in the 
overlapping of stability regions of the L = 1 and L = 2 phases^. Also near the other transition points the states with 
P > 1 could become stable so transitions could split into series of first order ones with long-periodic intermediate 
phases. Yet, all this cannot change the main conclusion about the infinite series of transitions accompanying the 
growth of interlayer disorder. 
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Note also that in the immediate vicinity of transition point T = J the present model has temperature dependencies 
of homogeneous magnetization and susceptibility characteristic to the ordinary mean- field ferromagnet. This is the 
consequence of the absence of disorder inside the layers undergoing purely ferromagnetic transitions. 



VI. CONCLUSIONS 



The present results give at least qualitative picture of magnetic structure and thermodynamic properties of layered 
Ising magnets having random interlayer exchange. Unlike the specific case considered here, really there always is 
some nonzero average interlayer interaction. It would split somewhat the potentials of degenerate glassy states, yet 
the growth of interlayer exchange fluctuations could still cause the changes of magnetic structure similar to ones 
described here. The study of such more realistic case can be done on the basis of present considerations. Qualitatively 
the presence of multiple almost degenerate metastable states with different total magnetizations would cause the 
inergodicity manifesting itself in a difference of 'field-cooled' and 'zero- field-cooled' parameters, specific hysteresis 
loop forms and long-time magnetic relaxation, which are really observed in many layered magnetsiiSi^i^i^. 

It should be also noted that the existence of a number of glassy phases at different u values means that transitions 
between them can also take place under temperature variations as u depends on T . Probably the observed in 
temperature oscillations of magnetization are the signatures of such transitions. One may also suppose that analogous 
yet more complex situation exists in case of short-range spin glasses, where a continuous changes of magnetic structure 
under temperature variations has been supposediS. 

One may also consider the present model as the intermediate variant between the mean-field Sherrington - Kirk- 
patrick spin glass and the short-range Edwards- Anderson model^. Combining the most of their essential features 
it allows for the first time to establish explicitly the existence of a number of metastable states and to study their 
properties which were not achieved in these classic models - there are only the estimates of a number of metastable 
states in the Sherrington - Kirkpatrick modeUi and their phenomenological description in the short-range spin glassi^. 
Now it is hard to say to what extent the existence of metastable states relies on the presence of long-range interac- 
tions. One may only assume, taking into account the evidences of their appearance in real disordered magnets with 
dipole-dipole and magnetoelastic interactions, that sufficiently slow decaying interactions ensures the existence of nu- 
merous metastable states. Anyway, the spin glass models without them have little chances to describe the inergodic 
phenomena observed in experiments. 
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